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INTRODUCTION 


For  over  seventy  years,  the  most  effective  and  widely  used  systemic  therapy  for  prostate  cancer 
has  been  androgen  deprivation  therapy  (ADT).  Although  nearly  all  men  will  respond  to  ADT, 
most  will  develop  castration-resistant  prostate  cancer  (CRPC).  However,  the  duration  of 
response  to  ADT  varies  widely.  Although  numerous  resistance  mechanisms  have  been 
described,  knowledge  of  them  typically  cannot  be  used  to  predict  a  priori  how  well  a  man  is 
likely  to  respond  to  ADT,  and  therefore  cannot  be  used  to  guide  management.  In  2013  we 
described  a  novel  resistance  mechanism  resulting  from  a  mutation  in  HSD3B1,  the  gene 
encoding  3P-hydroxysteroid  dehydrogenase- 1  (3f]HSD  1 ),  the  enzyme  responsible  for  catalyzing 
the  rate-limiting  step  in  the  conversion  of  adrenal  androgen  precursors  to  dihydrotestosterone 
(DHT),  and  which  is  required  for  all  pathways  of  DHT  synthesis.  HSD3B1  (1245 A>C)  changes 
amino  acid  367N— >T  and  renders  3f:SHSD  I  resistant  to  proteasomal  degradation,  causing 
profound  accumulation  and,  effectively,  gain-of-function.  The  resultant  increased  intratumoral 
metabolic  flux  of  adrenal  precursors  to  more  potent  androgens,  including  DHT  (the  most  potent 
androgen),  thus  enhances  androgen  receptor  activation  and  accelerates  tumor  proliferation, 
despite  castration.  Notably,  while  the  HSD3B i(1245C)  allele  can  be  acquired  by  mutation,  it  is 
also  heritable  in  the  form  of  a  single  nucleotide  polymorphism.  Recognizing  the  possibility  for  a 
novel  prostate  cancer  biomarker  with  the  potential  to  predict  resistance  to  ADT,  we  hypothesized 
that  men  inheriting  the  HSD3B1(1245C)  allele  would  exhibit  resistance  to  ADT.  Our  preliminary 
data  in  a  cohort  of  men  with  biochemical  failure  strongly  supported  this  hypothesis,  and 
prompted  the  current  project  to  validate  and  expand  these  findings. 
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ACCOMPLISHMENTS 


What  are  the  major  goals  of  the  project?  What  was  accomplished  under  these  goals? 

Major  goals  include  both  research  and  training  objectives.  The  research  objectives  are  discussed 
here,  and  the  training  goals  are  described  in  the  item  below  on  training  opportunities  and 
professional  development. 

Specific  Aim  1 :  validate  the  predictive  value  of  HSD3B1  genotype  in  independent  cohorts.  We 
have  worked  extensively  on  this  goal,  and  recently  published  this  work  in  Lancet  Oncology.  Our 
target  time  frame  was  14  months  for  the  first  milestone:  validation  of  the  predictive  value  of 
HSD3B1  genotype  in  an  independent  cohort  using  the  Mayo  Clinic  data.  We  completed  this  task 
in  12  months. 

Subtask  1:  Define  the  Mayo  Clinic  validation  cohorts.  The  post-prostatectomy  validation 
cohort  was  obtained  from  the  Sponsored  Programs  of  Research  Excellence  (SPORE) 
database  at  Mayo  Clinic.  We  identified  all  men  who:  (1)  received  ADT  for  biochemical 
and/or  non-metastatic  clinical  failure  following  prostatectomy,  and  (2)  had  blood  samples 
available  for  research.  To  ascertain  the  impact  of  HSD3B1  genotype  among  men  with  more 
advanced  disease,  we  queried  a  separate  Mayo  database  and  identified  men  enrolled  with 
metastatic  CRPC  who  had  blood  samples  banked  for  research.  To  minimize  confounding, 
patients  were  excluded  if  they  had  received  radiation  therapy  for  initial  local  therapy,  since 
use  of  ADT  with  RT  was  not  adequately  captured  in  this  cohort.  Patients  with  confirmed  or 
probable  neuroendocrine  prostate  cancer  were  excluded,  given  that  such  cancers  are 
generally  not  dependent  on  the  AR  axis. 

Subtask  2:  Obtain  germline  DNA  for  each  of  the  men  in  the  validation  cohorts  and  then 
determine  HSD3B1  genotype  using  a  high-resolution  melting  developed  for  this  purpose. 
We  obtained  peripheral  blood  mononuclear  cells  for  genotyping  in  the  Mayo  cohorts. 
Notably,  HSD3B1  has  one  homolog  ( HSD3B2 )  and  four  non-processed  pseudogenes,  which 
have  very  closely  related  DNA  sequences  that  may  obscure  detection  of  the  variant 
sequence.  We  therefore  developed  a  melting  assay  using  an  unlabeled  locked  nucleic  acid 
oligonucleotide  probe  in  an  asymmetric  polymerase  chain  reaction  to  perform  targeted 
genotyping  at  this  locus. 

DNA  from  the  Mayo  cohorts  was  extracted  from  peripheral  blood  mononuclear  cells. 
The  DNA  quantity  and  quality  was  measured  using  NanoDrop  (Thermo  Scientific, 
Waltham,  MA,  USA).  Genomic  DNA  was  stored  at  no  more  than  4°C.  The  PCR  was 
performed  in  20j.iL  reactions  (final  volume)  using  LightScanner  Master  Mix  (Biofire,  Salt 
Lake  City,  UT,  USA),  4pL  extracted  genomic  DNA  ( 1  Oug/pL),  forward  primer,  reverse 
primer,  and  unlabeled  LNA  probe  at  concentrations  of  0  05,  0-5,  and  0-5pM,  respectively 
(Table  SI  in  the  appendix  to  the  manuscript).  The  PCR  and  the  melting  analyses  were 
performed  in  StepOnePlus™  System  (Life  Technologies,  Grand  Island,  NY,  USA)  using  the 
following  conditions:  2  min  at  95°C,  followed  by  60  cycles  of  94°C  for  30  s,  66°C  for  30  s, 
and  75°C  for  30  s.  The  amplification  cycles  were  then  followed  by  the  melting  steps:  95°C 
for  30  s,  followed  by  cooling  to  25°C  for  60  s,  and  a  slow  final  denaturation  to  95°C  at  a 
thennal  ramp  rate  of  0  3%.  Data  were  analyzed  using  probe  melting  derivative  peaks.  Each 
validation  run  was  performed  with  duplicate  control  samples  representing  the  three  different 
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genotypes  and  each  sample  was  tested  in  duplicate  during  the  study  genotyping.  Our 
melting  analysis  technique  was  able  to  reliably  distinguish  genotypes  (see  Figure  SI  in  the 
appendix  to  the  manuscript),  and  was  validated  with  Sanger  sequencing  for  60  samples  (20 
for  each  genotype)  with  100%  concordance. 

Subtask  3:  Extract  demographic,  treatment,  and  clinical  outcomes  data  from  the  clinical 
annotation  in  Mayo  Clinic’s  prostate  cancer  registries.  For  demographic  and  treatment 
factors,  please  see  Table  1  in  the  manuscript  and  Tables  S2  and  S3  in  the  appendix  to  the 
manuscript.  Note  that  Table  1  also  includes  the  demographic  and  treatment  factors  of  the 
Cleveland  Clinic  primary  cohort,  which  was  the  cohort  used  for  the  preliminary  data  guiding 
this  project.  Outcomes  data  on  clinical  events  were  likewise  extracted — see  Sub  task  4  for 
additional  information  on  outcomes. 

Subtask  4:  Statistical  analysis.  Progression-free  survival,  distant  metastasis-free  survival, 
and  overall  survival  were  analyzed  using  Kaplan-Meier  methods.  Log-rank  tests  were  used 
to  compare  genotypes.  Cox  proportional  hazards  regression  was  perfonned  to  further  assess 
potential  allele-dosage  effects.  Demographic  and  treatment  characteristics  were  compared 
across  genotypes  to  assess  for  confounders  using  Fisher’s  exact  test  and  Kruskal-Wallis 
analysis  of  variance.  Fisher’s  exact  test  was  utilized  rather  than  the  Chi-square  test  because 
the  latter  is  less  accurate  when  counts  are  low.  Multivariable  Cox  proportional  hazards 
models  were  used  to  account  for  potential  confounders.  All  tests  were  two-sided,  and  P- 
values  less  than  or  equal  to  0-05  were  interpreted  as  statistically  significant.  Analyses  were 
performed  with  the  use  of  SAS  software,  version  9-3  (SAS  Institute,  Cary,  NC).  Time-to- 
event  outcomes  were  measured  from  the  date  ADT  was  originally  initiated.  The  key  results 
are  shown  in  Figures  1-3  in  the  manuscript,  which  demonstrate  that  progression- free 
survival,  distant  metastasis-free  survival,  and  overall  survival  varied  according  to  HSD3B1 
genotype.  More  specifically,  outcomes  worsened  as  the  number  of  inherited  variant  alleles 
increased.  Additional  details  are  included  in  the  text  of  the  Results  section  of  the 
manuscript  and  in  Tables  S2-S4  in  the  appendix  to  the  manuscript. 

Specific  Aim  2:  Perform  an  exploratory  analysis  of  serum  steroid  profiles  (delta-4/delta-5  ratios) 
to  evaluate  for  a  biochemical  correlate  to  HSD3B1  genotype,  which  could  represent  an  alternate 
biomarker.  Status:  to  be  done  (validation  was  the  first  priority). 

Specific  Aim  3:  Detennine  whether  HSD3B1  genotype  retains  predictive  value  when  ADT  is 
combined  with  radiotherapy  for  prostate  cancer.  Status:  to  be  done  (validation  was  the  first 
priority). 

What  opportunities  for  training  and  professional  development  has  the  project  provided? 

One  of  the  major  goals  of  this  project  is  to  foster  my  training  and  educational  development  in 
prostate  cancer  research,  and  in  this  regard  I  have  been  very  fortunate.  Now  at  a  critical  point  in 
my  development  as  a  young  investigator,  Dr.  Nima  Sharifi’s  ongoing  mentorship  has  been 
invaluable.  I  have  learned  much  through  his  modeling  of  the  scientific,  leadership,  professional, 
and  practical  elements  necessary  for  an  investigator  to  contribute  to  the  field.  I  have  also 
benefited  from  the  mentorship  of  Dr.  Felix  Feng,  who  has  likewise  been  a  role  model  and  source 
of  valuable  advice. 
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Apart  from  mentorship,  I  have  benefited  from  ongoing  educational  development  and 
training  through  my  involvement  in  multidisciplinary  genitourinary  conferences,  prostate  cancer 
research  seminars,  visiting  professor  sessions,  radiation  oncology  grand  rounds,  and  meetings 
with  Dr.  Sharifi  and  collaborators.  Additionally,  I  have  had  the  privilege  of  presenting  our  work 
at  the  2016  ASCO  annual  meeting  and  the  2016  Department  of  Defense  IMPaCT  meeting.  I  also 
look  forward  to  presenting  our  results  at  the  2016  Prostate  Cancer  Foundation  scientific  retreat. 
These  meeting  have  furnished  excellent  opportunities  to  network  with  other  prostate  cancer 
investigators  and  stay  abreast  of  developments  in  the  field. 

Finally,  I  have  had  the  privilege  of  caring  for  my  own  patients  with  prostate  cancer, 
which  has  yielded  experience  and  perspective  that  cannot  be  gained  any  other  way,  and  which 
continually  reaffirms  my  pledge  to  improve  care  for  future  men  with  prostate  cancer. 

How  were  the  results  disseminated  to  communities  of  interest? 

Please  see  above  regarding  presentation  at  national  meetings.  Please  see  Appendices  1  and  2  for 
published  results. 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

We  plan  to  work  on  Specific  Aim  2  listed  above,  in  an  effort  evaluate  whether  serum  steroid 
profiles  (delta-4/delta-5  ratios)  correlate  with  HSD3B1  genotype  and  whether  such  profiles  could 
represent  an  alternative  biomarker. 
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IMPACT 


What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Please  see  “Research  in  context”  box  on  the  second  page  of  the  journal  article  in  Appendix  1.  In 
brief,  we  assessed  whether  there  is  a  clinical  impact  of  inheriting  a  variant  form  (allele)  of  the 
HSD3B1  gene,  previously  shown  to  enable  prostate  cancer  cells  to  make  their  own  derivative  of 
testosterone  that  can  drive  tumor  progression,  despite  castration.  We  have  demonstrated  in  three 
independent  groups  of  men  treated  with  androgen  deprivation  therapy  (ADT)  that  the  number  of 
inherited  copies  of  the  variant  allele  corresponds  to  clinical  outcomes  (i.e.  likelihood  of  disease 
progression,  risk  of  spread  to  the  bones  and  other  distant  sites,  and  death).  Men  possessing  the 
variant  HSD3B1  allele  have  worse  clinical  outcomes  than  men  who  have  not  inherited  that  allele. 
Moreover,  men  who  have  two  copies  of  the  variant  allele  tend  to  have  profound  resistance  to 
ADT.  These  results  indicate  that  HSD3B1  genotype  is  a  powerful  genetic  biomarker  of 
resistance  to  ADT,  and  can  be  used  to  identify  men  who  might  benefit  from  escalated  therapy. 
Identification  of  these  men  has  potential  ramifications  for  guiding  clinicians  in  augmenting 
standard  ADT  with  highly-potent  inhibitors  of  the  androgen  receptor  pathway  or  chemotherapy 
in  an  effort  to  improve  outcomes.  Finally,  design  and  analysis  of  future  studies  should  take  into 
account  HSD3B1  genotype,  in  view  of  the  differences  in  outcomes  according  to  the  number  of 
variant  alleles  a  man  inherits. 

What  was  the  impact  on  other  disciplines? 

Nothing  to  report. 

What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report. 


CHANGES/PROBLEMS 

Nothing  to  report. 
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HSD3B1  and  resistance  to  androgen-deprivation  therapy  in 
prostate  cancer:  a  retrospective,  multicohort  study 

Jason  W  D  Hearn*,  Ghada  AbuAli*,  Chad  A  R eichard,  Chandana  A  Reddy,  Cristina  Magi-Galluzzi,  Kai-Hsiung  Chang,  Rachel  Carlson, 

Laureano  Rangel,  Kevin  Reagan,  Brian  J  Davis,  R  Jeffrey  Karnes,  Manish  Kohli,  Donald  Tindall,  Eric  A  Klein,  NimaSharifi 

Summary 

Background  HSD3B1  (1245A>C)  has  been  mechanistically  linked  to  castration-resistant  prostate  cancer  because  it 
encodes  an  altered  enzyme  that  augments  dihydrotestosterone  synthesis  from  non-gonadal  precursors.  We  postulated 
that  men  inheriting  the  HSD3B1  (1245C)  allele  would  exhibit  resistance  to  androgen-deprivation  therapy  (ADT). 

Methods  In  this  multicohort  study,  we  determined  HSD3B1  genotype  retrospectively  in  men  treated  with  ADT  for 
post-prostatectomy  biochemical  failure  and  correlated  genotype  with  long-term  clinical  outcomes.  We  used  data  and 
samples  from  prospectively  maintained  prostate  cancer  registries  at  the  Cleveland  Clinic  (Cleveland,  OH,  USA; 
primary  study  cohort)  and  the  Mayo  Clinic  (Rochester,  MN,  USA;  post-prostatectomy  and  metastatic  validation 
cohorts).  In  the  post-prostatectomy  cohorts,  patients  of  any  age  were  eligible  if  they  underwent  prostatectomy  between 
Jan  1, 1996,  and  Dec  31,  2009  (at  the  Cleveland  Clinic;  primary  cohort),  or  between  Jan  1, 1987,  and  Dec  31,  2011  (at  the 
Mayo  Clinic;  post-prostatectomy  cohort)  and  were  treated  with  ADT  for  biochemical  failure  or  for  non-metastatic 
clinical  failure.  In  the  metastatic  validation  cohort,  patients  of  any  age  were  eligible  if  they  were  enrolled  at  Mayo 
Clinic  between  Sept  1,  2009,  and  July  31,  2013,  with  metastatic  castration-resistant  prostate  cancer.  The  primary 
endpoint  was  progression-free  survival  according  to  HSD3B1  genotype.  We  did  prespecified  multivariable  analyses  to 
assess  the  independent  predictive  value  of  HSD3B1  genotype  on  outcomes. 

Findings  We  included  and  genotyped  443  patients:  118  in  the  primary  cohort  (who  underwent  prostatectomy),  137  in 
the  post-prostatectomy  validation  cohort,  and  188  in  the  metastatic  validation  cohort.  In  the  primary  study  cohort, 
median  progression-free  survival  diminished  as  a  function  of  the  number  of  variant  alleles  inherited:  6-6  years 
(95%  Cl  3  ■  8-not  reached)  in  men  with  homozygous  wild-type  genotype,  4  ■  1  years  (3  •  0-5  •  5)  in  men  with  heterozygous 
variant  genotype,  and  2  ■  5  years  (0  •  7  to  not  reached)  in  men  with  homozygous  variant  genotype  (p=0  ■  Oil) .  Relative  to 
the  homozygous  wild-type  genotype,  inheritance  of  two  copies  of  the  variant  allele  was  predictive  of  decreased 
progression-free  survival  (hazard  ratio  [HR]  2  -4  [95%  Cl  11-5  -3],  p=0-029),  as  was  inheritance  of  one  copy  of  the 
variant  allele  (HR  1-7  [1-0-2  -  9],  p=0-041).  Findings  were  similar  for  distant  metastasis-free  survival  and  overall 
survival.  The  effect  of  the  HSD3B1  genotype  was  independently  confirmed  in  the  validation  cohorts. 


Interpretation  Inheritance  of  the  HSD3B1  (1245C)  allele  that  enhances  dihydrotestosterone  synthesis  is  associated 
with  prostate  cancer  resistance  to  ADT.  HSD3B1  could  therefore  potentially  be  a  powerful  genetic  biomarker  capable 
of  distinguishing  men  who  are  a  priori  likely  to  fare  favourably  with  ADT  from  those  who  harbour  disease  liable  to 
behave  more  aggressively,  and  who  therefore  might  warrant  early  escalated  therapy. 


Funding  Prostate  Cancer  Foundation,  National  Institutes  of  Health,  US  Department  of  Defense,  Howard  Hughes 
Medical  Institute,  American  Cancer  Society,  Conquer  Cancer  Foundation  of  the  American  Society  of  Clinical 
Oncology,  Cleveland  Clinic  Research  Programs  Committee  and  Department  of  Radiation  Oncology,  Gail  and  Joseph 
Gassner  Development  Funds. 


Introduction 

Nearly  all  prostate  cancers  express  the  androgen  receptor, 
the  importance  of  which  is  underscored  by  androgen- 
deprivation  therapy  (ADT),  the  most  effective  and  widely 
used  systemic  therapy  for  prostate  cancer  for  the  past 
70  years.1  When  combined  with  radiotherapy,  ADT 
improves  survival  in  selected  patients.2,3  Similarly,  ADT 
confers  a  survival  advantage  if  given  immediately  after 
prostatectomy  in  node-positive  disease.4  ADT  is  the 
cornerstone  of  treatment  in  men  with  metastatic  disease,5 
and  has  shown  benefit  even  in  the  setting  of  biochemical 
failure  after  local  therapy  (ie,  in  the  context  of  increasing 
prostate-specific  antigen  [PSA]  after  definitive  surgery  or 


radiotherapy).6  For  example,  the  TOAD  trial7  published  in 
2016  showed  improved  survival  with  early  versus  delayed 
ADT  in  non-metastatic  men  with  increasing  PSA  level, 
most  of  whom  had  biochemical  failure  after  local  therapy. 
Although  nearly  all  men  initially  show  a  response  to 
ADT,  most  will  eventually  develop  castration-resistant 
prostate  cancer.  However,  the  duration  of  response  to 
ADT  varies  substantially.8,9 

Progression  from  castration-sensitive  to  castration- 
resistant  prostate  cancer  hinges  on  androgen  receptor 
reactivation,  which  can  occur  by  several  mechanisms.10,11 
Since  the  underlying  processes  typically  emerge  under 
selection  pressure  from  ADT,  knowledge  of  such  processes 
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Research  in  context 
Evidence  before  this  study 

We  searched  PubMed  for  studies  published  between 
Jan  1, 1980,  and  Feb  1,  2016,  on  biomarkers  of  resistance  to 
androgen-deprivation  therapy  (ADT)  in  prostate  cancer. 

We  used  the  search  terms  "androgen  deprivation  therapy" 

AND  (“resistance"  or  "efficacy"),  AND  "prostate  cancer",  AND 
“biomarker".  We  identified  144  studies,  common  limitations 
of  which  included  dependence  on  information  not  available  at 
ADT  initiation,  emphasis  on  associations  between 
polymorphisms  with  no  known  mechanistic  correlate,  absence 
of  association  with  endpoints  other  than  time  to  progression, 
and  absence  of  validation.  We  restricted  our  search  to  studies 
that  reported  on  people  and  were  published  in  English. 

Added  value  of  this  study 

We  tested  the  prespecified,  mechanism-based  hypothesis  that 
men  inheriting  the  HSD3B1  (1245C)  allele  would  exhibit 
innate  resistance  to  ADT,  since  this  allele  had  been  established 
to  increase  intratumoural  conversion  of  androgen  precursor 
steroids  to  more  potent  androgens  that  can  drive  disease 
progression,  despite  castration.  We  demonstrated  that 
HSD3B1  genotype  has  a  large  effect  on  long-term  clinical 


endpoints,  including  progression-free  survival,  distant 
metastasis-free  survival,  and  overall  survival.  Moreover,  a 
stepwise  difference  was  noted  in  outcomes  according  to  the 
number  of  HSD3B1  (12450  alleles  inherited,  with  differences 
measured  in  years.  Additionally,  we  have  validated  the 
influence  of  HSD3B1  genotype  for  progression-free  survival  by 
analysing  three  independent  cohorts,  representing  both  the 
post-prostatectomy  and  metastatic  disease  contexts.  Finally, 
we  found  that  HSD3B1  genotype  was  associated  with  overall 
survival  in  both  study  cohorts  (primary  and  metastatic 
validation  cohorts)  with  mature  follow-up. 

Implications  of  all  the  available  evidence 

Our  data  show  that  HSD3B1  genotype  is  a  powerful  genetic 
biomarker  of  resistance  to  ADT,  and  can  be  used  to  identify 
men  who  might  benefit  a  priori  from  escalated  therapy. 

This  identification  might  have  ramifications  for  selective  early 
incorporation  of  highly-potent  inhibitors  of  the  androgen 
receptor  axis,  and  could  potentially  inform  selection  of 
patients  for  chemohormonal  therapy.  Future  studies  should 
stratify  by  HSD3B1  genotype  in  view  of  the  differences  in 
outcomes  according  to  the  number  of  variant  alleles  present. 


generally  cannot  be  used  to  establish  a  priori  how  patients 
will  respond  to  this  treatment.  A  major  advance  in  the  past 
decade  has  been  the  increased  appreciation  of  intra¬ 
tumoural  androgen  synthesis.12'13  Before  ADT,  tumour 
androgen  supply  is  dominated  by  gonadal  testosterone. 
With  gonadal  suppression  during  ADT,  the  serum 
concentration  of  testosterone  is  greatly  depleted, 
inhibiting  tumour  growth.  However,  proliferation  can 
continue  in  the  context  of  intratumoural  androgen 
synthesis,  mostly  from  adrenal  precursor  steroids  and, 
possibly,  at  least  in  part,  due  to  de-novo  synthesis  from 
cholesterol.  Strong  evidence  for  the  importance  of 
intratumoural  androgen  synthesis  is  found  in  the  survival 
benefit  from  abiraterone  (which  depletes  intratumoural 
androgens)  and  enzalutamide  (which  competes  with 
intratumoural  androgens). 14-18  Additionally,  transcripts  for 
several  steroidogenic  enzymes — including  AKR1C3, 
HSD3B1,  and  HSD3B2 — are  consistently  upregulated  in 
castration-resistant  prostate  cancer.19  In  2013,  a  mutation 
in  HSD3B1  was  shown  to  provide  a  novel  mechanism  of 
resistance  to  ADT.20  HSD3B1  encodes  3fS-hydroxysteroid 
dehydrogenase-1  (3|3HSD1),  an  isoenzyme  that  is  mainly 
expressed  in  peripheral  tissues  (eg,  the  prostate,  skin, 
breast,  and  placenta),  is  responsible  for  catalysing  the 
rate-limiting  step  in  the  conversion  of  adrenal  androgen 
precursors  to  dihydrotestosterone,  and  is  required  for  all 
pathways  of  dihydrotestosterone  synthesis.21  HSD3B1 
(1245A>C)  changes  aminoacid  367Asn-»Thr  and  renders 
3(3HSD1  resistant  to  proteasomal  degradation,  causing 
substantial  accumulation  of  this  enzyme  and,  effectively, 
gain-of-function.  The  resultant  increased  intratumoural 
metabolic  flux  of  adrenal  precursors  to  more  potent 


androgens,  such  as  dihydrotestosterone  (the  most  potent 
androgen),  therefore  enhances  androgen  receptor 
activation  and  accelerates  tumour  proliferation,  despite 
castration.  Although  the  HSD3B1  (1245C)  allele  can  be 
acquired  by  mutation,  it  is  also  heritable  in  the  form  of  a 
single  nucleotide  polymorphism  (SNP;  rsl047303)  with 
allelic  frequency  of  15-35%  in  most  cohorts  (this  frequency 
varies  by  ethnic  origin,  and  is  much  lower  in  Asians  and 
higher  in  white  people).22 

At  present,  the  clinical  relevance  of  HSD3B1  (1245C) 
inheritance  in  prostate  cancer  is  unknown.  One  study23 
in  China  has  linked  this  variant  allele  with  increased 
progression  to  castration-resistant  prostate  cancer,  but 
not  with  other  endpoints.  However,  that  study23  was 
limited  by  the  extremely  low  prevalence  of  the  HSD3B1 
(1245C)  allele  in  Chinese  men,  and  therefore  included  no 
men  homozygous  for  the  variant  allele  and  very  few 
who  had  even  one  allele.  Additionally,  the  study  was 
confounded  by  the  fact  that  nearly  twice  as  many  men 
included  in  the  homozygous  wild-type  cohort  started 
ADT  when  they  already  had  metastatic  disease.23  As  such, 
the  potential  clinical  relevance  of  the  variant  allele 
remained  undetermined. 

Our  present  study  was  designed  to  specifically  test  the 
mechanism-based  hypothesis  that  men  inheriting 
the  variant  HSD3B1  (1245C)  allele  would  demonstrate 
evidence  of  intrinsic  resistance  to  ADT,  which  would 
manifest  as  shorter  time  to  progression,  worse  distant 
metastasis-free  survival,  and  potentially  worse  overall 
survival  than  men  not  inheriting  this  variant.  Moreover, 
on  the  basis  of  the  number  of  variant  alleles  inherited, 
we  hypothesised  that  men  homozygous  for  the  variant 
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allele  would  have  the  worst  outcomes,  men  homozygous 
for  the  wild-type  allele  would  have  the  best  outcomes, 
and  heterozygotes  would  potentially  have  intermediate 
outcomes. 

Methods 

Study  design  and  participants 

In  this  multicohort  study,  we  used  data  from  three  large, 
prospectively  maintained  prostate  cancer  registries:  the 
primary  study  cohort  at  Cleveland  Clinic  (Cleveland,  OH, 
USA)  and  two  validation  cohorts  from  the  Mayo  Clinic 
(Rochester,  MN,  USA).  For  eligibility  in  the  primary  study 
cohort,  men  (of  any  age)  had  to  have  undergone 
prostatectomy  at  the  Cleveland  Clinic  between  Jan  1,  and 
Dec  31,  2009,  and  to  have  been  treated  with  ADT  for 
biochemical  failure.  Biochemical  failure  was  defined  as  a 
PSA  value  of  0-2  pg/L  or  higher,  followed  by  a  second 
higher  reading  (no  specified  time  interval  between  the 
two  readings),  or  one  PSA  value  of  0-5  pg/L  or  higher.24 
We  specified  this  definition  to  reduce  the  likelihood 
of  inadvertently  including  men  with  low  detectable 
concentration  of  PSA  produced  from  residual  benign 
tissue  or  biologically  indolent  cancer,  since  the  inclusion  of 
such  patients  would  undermine  the  study’s  ability  to  test 
the  mechanism-based  hypothesis.  Patients  with  minimally 
detectable  PSA  coming  from  benign  tissue  rather  than 
from  cancer  would  not  have  progression,  metastasis,  or 
death  from  cancer  (because  they  had  no  cancer).  Patients 
with  a  minute  amount  of  very  indolent  cancer  such  that 
they  never  met  the  definition  of  biochemical  failure  would 
be  less  likely  to  progress,  metastasise,  or  die  from  cancer 
than  those  whose  cancer  proliferated  fast  enough  that  they 
met  the  definition  of  biochemical  failure. 

Additionally,  men  who  received  postoperative  adjuvant 
or  salvage  radiotherapy  were  eligible  for  the  study,  if  they 
had  residual  active  disease  as  indicated  by  continued 
increase  in  their  PSA  after  treatment  (ie,  PSA  after 
radiotherapy  higher  than  postoperative  preradiotherapy 
value).  Patients  were  excluded  if  the  quality  of  their 
sample  of  DNA  was  insufficient  to  permit  genotyping. 

For  the  post-prostatectomy  validation  cohort,  we  used 
patient  data  from  the  National  Cancer  Institute-funded 
Sponsored  Programs  of  Research  Excellence’s  (SPORE) 
database  at  the  Mayo  Clinic,  and  included  men  who 
underwent  prostatectomy  between  Jan  1,  1987,  and 
Dec  31,  2011.  We  identified  and  included  all  men  who 
had  received  ADT  for  biochemical  or  non-metastatic 
clinical  failure  (eg,  local  failure)  after  prostatectomy  and 
had  blood  samples  available  for  research.  For  the 
metastatic  validation  cohort,  we  queried  a  large  Mayo 
Clinic  registry  of  patients  with  advanced  prostate  cancer. 
We  identified  men  enrolled  between  Sept  1,  2009,  and 
July  31,  2013,  with  metastatic  castration-resistant  prostate 
cancer  who  had  blood  samples  banked  for  research. 

To  minimise  confounding,  we  excluded  patients  from 
the  metastatic  validation  cohort  if  they  had  received 
radiotherapy  for  initial  local  therapy  (since  use  of  ADT 


with  radiotherapy  was  not  adequately  captured  in 
this  cohort)  or  if  they  had  confirmed  or  probable 
neuroendocrine  prostate  cancer  (given  that  such  cancers 
are  generally  not  dependent  on  the  androgen  receptor 
axis).  We  used  no  age  restrictions  for  patients  in  the 
validation  cohorts.  All  patients  were  aged  older  than 
18  years. 

All  biological  samples  and  clinical  data  were  obtained 
with  written  informed  patient  consent  as  part  of  informed 
consent  protocols  approved  by  local  institutional  review 
boards  at  the  Cleveland  Clinic  and  the  Mayo  Clinic. 

Procedures 

We  reviewed  prostatectomy  specimens  from  the  primary 
cohort  to  obtain  cores  of  non-neoplastic  tissue  for 
germline  genotyping  of  all  patients  in  the  primary  cohort. 

A  genitourinary  pathologist  (CM-G)  reviewed  each  of 
these  specimens  to  identify  non-neoplastic  tissue. 

Genotyping  was  accomplished  with  a  melting  assay 

(appendix  pp  1,  2,  7).  We  used  peripheral  blood  see  online  for  appendix 

mononuclear  cells  to  genotype  the  two  validation  cohorts 

using  the  melting  assay.  HSD3B1  has  one  homologue 

(HSD3B2)  and  four  non-processed  pseudogenes 

(HSD3BP1,  HSD3BP2,  HSD3BP4,  and  HSD3BP5),  which 

have  very  closely  related  DNA  sequences  that  might 

obscure  detection  of  the  variant  sequence.21  We  therefore 

developed  a  melting  assay  using  an  unlabelled,  locked, 

nucleic  acid  oligonucleotide  probe  in  an  asymmetric  PCR 

to  complete  targeted  genotyping  at  this  locus  for  all 

cohorts  (appendix  pp  1,  7).  All  patients  were  identified  in 

all  cohorts  (by  CARed,  JWDH,  RC,  LR,  and  MK)  before 

determination  of  genotype,  and  tissue  processing  and 

genotyping  were  completed  by  investigators  who  were 

masked  to  clinical  data.  Data  collection  and  analysis  were 

uncoupled,  because  the  study  statistician  (CARed)  was  not 

involved  in  collection  of  clinical  data  or  genotyping. 

Outcomes 

The  primary  endpoint  was  progression-free  survival. 

Progression  was  an  investigator-assessed  composite 
endpoint  defined  as  the  first  occurrence  of  a  second 
increase  in  PSA  on  ADT  (no  absolute  threshold), 
radiographic  or  clinical  progression  (defined  as  the 
development  of  a  local  or  regional  nodal  recurrence  on 
physical  examination,  CT,  or  MRI;  distant  metastasis 
seen  on  bone  scan,  CT,  or  MRI;  or  death  from  prostate 
cancer),  or  initiation  of  any  second-line  therapy.  The 
secondary  endpoints  were  distant  metastasis-free 
survival  (ie,  the  probability  of  surviving  without  distant 
metastases)  and  overall  survival.  Time-to-event  outcomes 
were  measured  from  the  date  when  ADT  was  initiated. 

Statistical  analysis 

We  analysed  progression-free  survival,  distant  metastasis- 
free  survival,  and  overall  survival  using  Kaplan-Meier 
methods.  We  used  the  log-rank  test  for  trend  with 
one  degree  of  freedom  to  assess  potential  gene-dose 
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effects  across  the  three  genotypes.  We  did  Cox 
proportional  hazards  regression  to  further  assess 
potential  allele-dose  effects.  We  compared  demographic 
and  treatment  characteristics  across  genotypes  to  assess 
for  confounders  using  Fisher’s  exact  test  and  the 
Kruskal-Wallis  one-way  analysis  of  variance.  We  used 
Fisher’s  exact  test  rather  than  the  x2  test:  because  the 
latter  is  less  accurate  if  the  number  of  patients  in  a 
category  is  low  (ie,  lower  than  five).  We  used  multivariable 
Cox  proportional  hazards  models  to  account  for  potential 


Homozygous  wild-type 
patients  (n=44) 

Heterozygous 
patients  (n=62) 

Homozygous  variant 
patients  (n=12) 

p  value* 

Follow-up,  years 

71  (40-90) 

6-7  (41-9-0) 

6-2  (3-6-6  8) 

0-36 

Age,  years 

65  (60-70) 

65  (55-70) 

63 (57-67) 

0-36 

Ethnic  origin 

0-16 

White 

34  (77%) 

56  (90%) 

11  (92%) 

Black 

8  (18%) 

4  (6%) 

0 

Hispanic 

0 

1(2%) 

0 

Asian 

1(2%) 

0 

0 

Other 

1(2%) 

1(2%) 

1  (8%) 

Pathological  T  stage! 

0-056 

T2 

14  (32%) 

9  (15%) 

0 

T3a 

13  (30%) 

17  (27%) 

4  (33%) 

T3b-4t 

17  (39%) 

36  (58%) 

8  (67%) 

Pathological  N  stage! 

0-97 

NO 

30  (68%) 

39  (63%) 

9  (75%) 

N1 

10  (23%) 

16  (26%) 

2  (17%) 

NX 

4  (9%) 

7  (11%) 

1  (8%) 

Prostate-specific  antigen 
at  ADT  initiation  (pg/L) 

3-3  (0-8-80) 

1-5  (0-6-54) 

1-9  (11-7-3) 

0-33 

Gleason  score 

0-42 

6 

2  (5%) 

0 

0 

7 

22  (50%) 

35  (57%) 

5  (42%) 

8-10 

20  (46%) 

27  (44%) 

7(58%) 

ADT  type 

0-69 

GnRH  agonist  or 
orchiectomy 

41  (93%) 

56  (90%) 

12  (100%) 

Androgen  receptor 
antagonist 

3  (7%) 

6  (10%) 

0 

ADT  use 

0-23 

Continuous 

25  (57%) 

27  (44%) 

8  (67%) 

Intermittent 

19  (43%) 

35  (57%) 

4(33%) 

Adjuvant  or  salvage  radiation  therapy 

0-66 

No 

31  (71%) 

38 (61%) 

8  (67%) 

Yes 

13  (30%) 

24  (39%) 

4(33%) 

Neoadjuvant  treatment 

0-97 

No 

31  (70%) 

45  (73%) 

9  (75%) 

ADT 

10  (23%) 

14  (23%) 

2  (17%) 

Chemotherapy  or 
immunotherapy 

3  (7%) 

3  (5%) 

1  (8%) 

Data  are  median  (IQR)  or  n  (%).  ADT=androgen-deprivation  therapy.  GnRH=gonadotropin-releasing  hormone. 

*p  values  are  comparisons  across  the  three  genotypes.  tPathological  stage  refers  to  the  American  Joint  Commission  on 
Cancer's  tumour-node-metastasis  staging  system.28  !StagesT3b  andT4  were  combined  because  only  one  patient 
(a  heterozygote)  was  stageT4  in  our  study  population. 

Table:  Baseline  characteristics  of  the  primary  cohort  by  HSD3B1  genotype 

confounders.  Covariates  included  common  clinico- 
pathological  variables  known  to  correlate  with  clinical 
outcomes,  including  American  Joint  Commission  on 
Cancer  T  stage,  N  stage,  PSA  at  initiation  of  ADT,  and 
Gleason  score.  We  examined  T  stage  as  a  categorical 
variable  using  all  T  stages  in  two  separate  models.  In  the 
first  model,  stage  T2  was  the  reference  group,  examining 
stage  3a  versus  stage  2  and  3b— 4-  versus  2.  In  the  second 
model,  stage  T3a  was  the  reference  group,  examining 
stage  3b — 4  versus  3a  and  2  versus  3a.  Stage  T4  was 
combined  with  stage  T3b  because  there  was  only 
one  stage  T4  patient  in  our  study  population.  Stage  T3b^I 
is  typically  more  predictive  of  poor  clinical  outcomes 
than  is  T3a,  and  therefore  we  did  not  combine  all 
T3  patients  into  one  category.  All  tests  were  two-sided, 
and  p  values  less  than  or  equal  to  0  •  05  were  judged  to  be 
statistically  significant.  We  used  SAS  version  9.3  for  our 
statistical  analyses. 

Role  of  the  funding  source 

The  funders  of  the  study  had  no  role  in  study  design, 
data  collection,  data  analysis,  data  interpretation,  or 
writing  of  the  report.  JWDH,  GA,  CARei,  CARed,  CM-G, 
RC,  LR,  MK,  DT,  and  NS  had  access  to  all  the  raw  data. 
After  extensive  quality  assurance  and  data  verification,  all 
authors  vouched  for  the  completeness  and  integrity  of 
the  data  and  statistical  analysis.  The  corresponding 
author  had  full  access  to  all  data  and  had  the  final 
responsibility  for  the  decision  to  submit  for  publication. 

Results 

In  the  primary  Cleveland  Clinic  cohort,  we  identified 
177  eligible  men  who  met  study  criteria  and  underwent 
prostatectomy  between  Jan  1,  1996,  and  Dec  31,  2009. 
Of  these,  118  (67%)  were  successfully  genotyped  and  could 
be  analysed,  which  is  indicative  of  the  well-established 
challenge  of  genotyping  with  formalin-fixed  paraffin- 
embedded  tissue.25-27  Our  melting  analysis  technique  was 
able  to  reliably  distinguish  genotypes  (appendix  p  2),  and 
was  validated  with  Sanger  sequencing  for  60  samples 
(20  for  each  genotype)  with  100%  concordance. 

Of  the  genotyped  men  in  the  primary  cohort,  HSD3B1 
(1245C)  allelic  frequency  was  36%  (86  of  236  alleles); 
44  (37%)  of  118  men  were  homozygous  wild-type,  62  (53%) 
were  heterozygous,  and  12  (10%)  were  homozygous 
variant.  No  significant  differences  were  noted  in 
demographic  or  treatment  characteristics  of  participants 
by  genotype,  although  the  crude  proportions  of  T3b— 4 
tumours  and  Gleason  score  8-10  disease  were  highest  in 
homozygous  variant  men  (table).  Multivariable  analyses 
were  used  to  adjust  accordingly  for  these  differences. 

Of  140  patients  in  the  post-prostatectomy  validation 
cohort  (the  Mayo  Clinic  SPORE  registry)  who  met  study 
criteria,  137  (98%)  were  successfully  genotyped  and  could 
be  analysed.  HSD3B1  (1245C)  allelic  frequency  was  26% 
(70  of  274  alleles),  and  the  genotype  distribution  was  as 
follows:  77  (56%)  of  137  men  were  homozygous  wild-type, 
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50  (37%)  were  heterozygous,  and  ten  (7%)  were 
homozygous  variant  (for  demographic  and  treatment 
characteristics  see  appendix  p  8).  Of  204  patients  who 
met  study  criteria  and  were  in  the  metastatic  validation 
cohort  (Mayo  Clinic  metastatic  prostate  cancer  registry) 
between  Sept  1,  2009,  and  July  31,  2013,  188  (92%) 
patients  were  successfully  genotyped  and  could  be 
analysed.  These  men  had  initially  been  diagnosed  with 
prostate  cancer  between  Jan  1,  1983,  and  July  31,  2012. 
HSD3B1  (1245C)  allelic  frequency  was  27%  (101  of 
376  alleles),  and  the  genotype  distribution  was  as  follows: 
98  (52%)  of  188  men  were  homozygous  wild-type, 
79  (42%)  were  heterozygous,  and  11  (6%)  were 
homozygous  variant  (for  demographic  and  treatment 
characteristics  see  appendix  p  9).  Pooling  data  from  all 
three  cohorts,  the  overall  HSD3B1  (1245C)  allelic 
frequency  was  29%  (257  of  886  alleles),  and  the  pooled 
genotype  distribution  was  as  follows:  219  (49%)  of 
443  men  were  homozygous  wild-type,  191  (43%)  were 
heterozygous,  and  33  (7%)  were  homozygous  variant. 

In  both  validation  cohorts,  the  demographic  and 
treatment  characteristics  of  the  patients  did  not  differ 
according  to  genotype  (appendix  pp  8-9),  with  the 
exception  of  N-stage  in  the  Mayo  Clinic  post-prosta- 
tectomy  SPORE  cohort.  The  distribution  of  N1  disease  in 
this  cohort  was  five  (50%)  of  ten  men  in  the  homozygous 
variant  group,  compared  with  ten  (13%)  of  77  men  in  the 
homozygous  wild-type  and  eight  (16%)  of  50  men  in  the 
heterozygous  group  (p=0-025).  This  discrepancy  was 
adjusted  for  during  multivariable  analysis.  Median 
follow-up  in  this  cohort  (4-9  years  [IQR  2 -6-7 -6])  was 
shorter  than  for  the  primary  cohort  (6  ■  7  years  [3  •  9-8  •  9]; 
table)  and  the  Mayo  metastatic  validation  cohort 
(6-0  years  [3 -1-9 -2]),  and  was  especially  short  in  the 
homozygous  variant  group  (3-2  years  [1  -  8 — 4  -  6]). 
Differences  in  the  duration  of  follow-up  between  the 
three  genotypes  were  taken  into  account  with  the  use  of 
Cox  proportional  hazards  regression.  In  view  of  the  low 
number  of  metastatic  and  death  events  during  the  short 
observation  period  (appendix  p  3),  the  post-prostatectomy 
validation  cohort  was  used  to  validate  the  effect  of 
HSD3B1  genotype  on  only  progression-free  survival. 
By  contrast,  the  median  follow-up  from  ADT  initiation 
for  the  metastatic  validation  cohort  was  sufficient  to 
enable  analysis  of  both  progression-free  survival  and 
overall  survival.  Distant  metastasis-free  survival  was  not 
calculated  in  this  cohort,  since  these  were  patients  who 
already  had  metastasis.  Chemotherapy  was  commonly 
used  for  castration-resistant  prostate  cancer,  and  its  use 
was  similar  across  the  three  genotypes  in  the  metastatic 
cohort  (57%  [56  of  98]  in  homozygous  wild-type  men, 
61%  [48  of  79]  in  heterozygotes,  and  55%  [six  of  11] 
in  homozygous-variant  men  [p=0  -  83];  appendix  p  9). 
Docetaxel  was  used  in  106  (96%)  of  the  110  patients  who 
received  chemotherapy  in  this  cohort. 

Progression-free  survival  was  significantly  associated 
with  HSD3B1  genotype  in  the  primary  study  cohort 


A  Primary  cohort 


Number  at  risk 
(number  censored) 


Homozygous  wild-type  44(3) 

28(7) 

15(5) 

7(6) 

1(1) 

Heterozygous  62  (5) 

36(4) 

17(4) 

4(3) 

1(1) 

Homozygous  variant  12  (0) 

5(0) 

3(3) 

0(0) 

0(0) 

B  Post-prostatectomy  validation  cohort 


Number  at  risk 
(number  censored) 

Homozygous  wild-type  77(1)  58(4)  43(5)  33(1) 

Heterozygous  50  (3)  39  (5)  28  (5)  15  (1) 

Homozygous  variant  10  (0)  4  (1)  2  (0)  0  (0) 

C  Metastatic  validation  cohort 


26(5) 

13(2) 

0(0) 


17(7) 

8(2) 

0(0) 


8(1) 

5(1) 

0(0) 


6(5) 

2(2) 

0(0) 


Number  at  risk 
(number  censored) 


Time  (years) 


Homozygous  wild-type  98(0) 

63  (0) 

46(0) 

32(0) 

26(0) 

20(0) 

15(0) 

10  (0) 

Heterozygous  79  (0) 

50(0) 

29(0) 

20(0) 

13(0) 

11(0) 

9(0) 

8(0) 

Homozygous  variant  11(0) 

4(0) 

2(0) 

1(0) 

1(0) 

1(0) 

0(0) 

0(0) 

Figure  1:  Progression-free  survival  according  to  HSD3B1  genotype 

Progression-free  survival  according  to  H5D3B1  genotype  in  (A)  the  primary  study  cohort,  (B)the  post-prostatectomy 
validation  cohort,  and  (C)  the  metastatic  validation  cohort. 
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Figure  2:  Distant  metastasis-free  survival  in  the  primary  study  cohort 


(figure  1A),  and  diminished  as  a  function  of  the  number 
of  HSD3B1  (1245C)  alleles  inherited:  median  6-6  years 
(95%  Cl  3  -8-not  reached)  in  homozygous  wild-type  men; 
4-1  years  (3 -0-5 -5)  in  heterozygous  men;  and  2-5  years 
(0-7-not  reached)  in  homozygous  variant  men  (p=0-011). 
Relative  to  the  homozygous  wild-type  genotype,  inheritance 
of  two  copies  of  the  variant  allele  was  predictive  of 
decreased  progression-free  survival  (hazard  ratio  [HR]  2-4 
[95%  Cl  1-1-5 -3],  p=0-029).  Additionally,  inheritance  of 
one  copy  of  the  variant  allele  was  also  associated  with 
decreased  progression-free  survival  (HR  1-7  [95%  Cl 
1-0-2  -9],  p=0-041).  On  multivariable  analysis  (appendix 
p  10),  the  HR  for  progression  was  1-6  (95%  Cl  1-0-2-7, 
p=0-074)  for  men  with  at  least  one  variant  allele,  which 
compared  favourably  with  Gleason  score  (for  Gleason 
score  8-10  vs  6-7:  HR  1-3  [95%  Cl  0- 8-2-0],  p=0-31). 
The  distribution  of  each  type  of  progression  event 
(first  occurrence  of  a  second  increase  in  PSA  on  ADT, 
radiographic  or  clinical  progression,  or  initiation  of 
second-line  therapy)  is  described  in  the  appendix  (p  11). 

In  the  post-prostatectomy  validation  cohort,  pro¬ 
gression-free  survival  was  associated  with  HSD3B1 
genotype  (figure  IB).  Median  progression-free  survival 
was  3  ■  3  years  (95%  Cl  1  ■  9^1  •  9)  in  homozygous  wild-type 
men,  2-8  years  (2 -1-5-1)  in  heterozygous  men,  and 
0-9  years  (0-2-2 -6)  in  homozygous  variant  men 
(p=0-0022).  Homozygous  variant  men  had  an  increased 
risk  of  progression  compared  with  homozygous  wild-type 
men  (HR  3-4  [95%  Cl  1- 6-7-0],  p=0-0013),  whereas 
progression-free  survival  did  not  dilfer  significantly 
between  heterozygous  men  and  the  wild-type  group 
(HR  1-0  [0 - 7—1  •  7] ,  p=0-85).  The  effect  of  homozygous 
variant  genotype  on  progression-free  survival  persisted 
with  prespecified  adjustment  for  lymph  node  status 
(lymph  node-negative  vs  lymph  node-positive;  HR  2-7 
[95%  Cl  1-2-5  -9],  p=0-013),  whereas  heterozygotes  did 
not  differ  significantly  from  homozygous  wild-type  men 
(1-0  [0-6-1 -6],  p=0-98). 


Median  progression-free  survival  in  the  metastatic 
validation  cohort  was  also  associated  with  HSD3B1 
genotype  (figure  1C).  As  in  the  other  two  cohorts,  the 
largest  difference  was  between  homozygous  wild-type 
(median  1-8  years  [95%  Cl  1-2-2 -5])  and  homozygous 
variant  men  (0-8  years  [0-3-1 -6],  p=0-024).  The 
corresponding  HR  for  progression  was  2-0  (95%  Cl 
1  ■  1-3  •  8,  p=0  •  027) .  Men  with  heterozygous  genotype  had 
an  intermediate  progression-free  survival  (1-4  years 
[95%  Cl  1-0-1 -8],  HR  for  progression  1-1  [0-8-1- 5], 
p=0  -  38). 

Distant  metastasis-free  survival  was  assessable  only  in 
the  primary  study  cohort  because  of  inadequate  length  of 
follow-up  in  the  post-prostatectomy  validation  cohort  and 
the  presence  of  metastatic  disease  at  ADT  initiation  in 
the  metastatic  cohort.  Distant  metastasis-free  survival 
was  therefore  not  validated  in  any  external  cohorts. 
Distant  metastasis-free  survival  was  significantly 
associated  with  HSD3B1  genotype  in  the  primary  study 
cohort  and  decreased  according  to  the  number  of  variant 
alleles  inherited:  median  9-1  years  (95%  Cl  7-4— not 
reached)  in  homozygous  wild-type  men;  6  •  8  years 
(4- 3-7-4)  in  heterozygotes;  and  3-6  years  (1-0-7 -3)  in 
homozygous  variant  men  (p=0-014;  figure  2).  Compared 
with  the  homozygous  wild-type  genotype,  inheritance  of 
two  copies  of  the  variant  allele  was  predictive  of  worse 
distant  metastasis-free  survival  (HR  2  •  7,  95%  Cl  1  •  2-6  •  2, 
p=0-022).  The  corresponding  HR  for  men  with  hetero¬ 
zygous  genotype  (ie,  one  copy  of  the  variant  allele)  was 
1-7  (95%  Cl  1-0-2-8,  p=0-074).  The  effect  of  both 
homozygous  variant  genotype  and  heterozygous 
genotype  on  distant  metastasis-free  survival  was 
confirmed  on  multivariable  analysis  (appendix  p  10). 

Overall  survival  was  significantly  associated  with 
HSD3B1  genotype  in  both  of  the  cohorts  with  sufficiently 
mature  follow  up  to  enable  its  analysis  (the  primary  and 
metastatic  validation  cohorts).  Figure  3A  shows  a 
stepwise  reduction  in  overall  survival  as  a  function  of  the 
number  of  variant  alleles  inherited  in  the  primary  cohort. 
The  5-year  point  estimates  of  overall  survival  were  82% 
(95%  Cl  69-94)  in  homozygous  wild-type  men,  74% 
(62-85)  in  heterozygous  men,  and  58%  (30-86)  in 
homozygous  variant  men;  and  10-year  point  estimates  of 
overall  survival  were  55%  (95%  Cl  35-75)  in  homozygous 
wild-type  men,  35%  (21-49)  in  heterozygous  men,  and 
0%  (0-0)  in  homozygous  variant  men  (p=0-0064). 
Median  overall  survival  in  the  primary  cohort  was 
16-7  years  (95%  Cl  7-9-16-7)  in  homozygous  wild-type 
men,  7-4  years  (6-4— 9-3)  in  heterozygous  men,  and 
7-3  years  (3 -8-8-1)  in  homozygous  variant  men 
(p=0-0064).  Relative  to  the  homozygous  wild-type 
genotype,  inheritance  of  two  copies  of  the  variant  allele 
was  predictive  of  worse  overall  survival  (HR  3  ■  3  [95%  Cl 
1-3-8 -3],  p=0-013).  For  men  with  the  heterozygous 
genotype,  the  inheritance  of  one  copy  of  the  variant  allele 
was  also  predictive  of  worse  overall  survival  (HR  2-0 
[1- 1-3  -7],  p=0-036).  Multivariable  analysis  confirmed  the 
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effect  of  HSD3B1  genotype  on  overall  survival 
(appendix  p  10). 

Similarly,  in  the  metastatic  validation  cohort,  median 
overall  survival  from  the  time  of  ADT  initiation  was 
reduced  according  to  the  number  of  variant  alleles 
present.  The  5-year  point  estimates  of  overall  survival 
were  72%  (95%  Cl  63-82)  in  homozygous  wild-type  men, 
62%  (51-73)  in  heterozygous  variant  men,  and  46% 
(16-75)  in  homozygous  variant  men;  and  10-year  point 
estimates  of  overall  survival  were  49%  (95%  Cl  38-60)  in 
homozygous  wild-type  men,  25%  (14—36)  in  heterozygous 
men,  and  16%  (0-39)  in  homozygous  variant  men 
(p=0-0042).  Median  overall  survival  in  the  metastatic 
validation  cohort  was  9-7  years  (95%  Cl  6- 7-12-1) 
in  homozygous  wild-type  men,  6-8  years  (5 -2-8-0)  in 
heterozygous  variant  men,  and  4-6  years  (1-6-7 -5) 
in  homozygous  variant  men  (p=0-0042;  figure  3B). 
Compared  with  men  with  homozygous  wild-type 
genotype,  the  HR  for  death  was  2-5  (95%  Cl  1- 2-5-0, 
p=0-013)  in  homozygous  variant  men  and  1-5  (1-  0-2-1, 
p=0-036)  in  heterozygotes.  Outcomes  for  each  cohort 
with  all  genotypes  combined  are  in  the  appendix  (pp  4—6). 
Prostate  cancer-specific  survival  for  the  primary  cohort 
(which  had  data  for  cause  of  death)  is  shown  in 
the  appendix  (p  6).  As  with  overall  survival,  prostate 
cancer-specific  survival  diminished  according  to  HSD3B1 
inheritance  between  the  three  genotypes  (p=0-029). 
Compared  with  the  homozygous  wild-type  reference 
group,  the  HR  for  prostate  cancer- specific  mortality  was 
3-1  (95%  Cl  1-1-8 -7;  p=0-027)  in  homozygous  variant 
men,  and  was  1-7  (0-8-3 -3;  p=0 - 17)  in  heterozygotes. 
Overall  survival  could  not  be  meaningfully  analysed  in 
the  SPORE  post-prostatectomy  cohort  due  to  inadequate 
length  of  follow-up. 

Discussion 

Our  results  show  that  inheritance  of  the  HSD3B1 
(1245C)  allele  (homozygous  variant  and  heterozygotes) 
was  associated  with  decreased  progression-free  survival, 
distant  metastasis-free  survival,  and  overall  survival, 
with  differences  measured  in  years.  These  results  are 
biologically  credible,  because  tumour  cells  carrying  the 
HSD3B1  (1245C)  allele  are  able  to  more  efficiently 
produce  their  own  dihydrotestosterone  than  those 
without  this  allele.  The  observed  stepwise  decrease  in 
all  three  endpoints  would  be  unlikely  to  occur  by 
chance,  and  strongly  suggests  allele  dose-dependence. 
Additionally,  validation  of  our  initial  findings  in 
two  independent  cohorts  for  progression-free  survival 
supports  the  effect  of  HSD3B1  (1245C)  inheritance  on 
resistance  to  ADT.  Genotype  might  also  be  prognostic 
irrespective  of  ADT,  although  the  mechanistic  under¬ 
pinnings  suggest  that  the  associated  growth  advantage 
would  be  most  pronounced  with  ADT. 

As  an  extension  of  our  preclinical  work  that  identified 
an  association  between  HSD3B1  (1245 C)  and  increased 
3|3HSD1  enzyme  concentrations,  we  previously  assessed 


Figure  3:  Overall  survival 

Overall  survival  according  to  HSD3B1  genotype  in  (A)  the  primary  study  cohort  and  (B)  the  metastatic  validation  cohort. 


tumours  in  men  with  castration-resistant  prostate 
cancer.20  In  homozygous  variant  tumours,  we  noted  a 
robust  increase  in  3(3HSD1  enzyme  concentrations, 
whereas  we  did  not  record  any  such  clear  increase  in 
heterozygous  variant  tumours.20  However,  expression 
level  was  measured  in  only  two  of  these  tumours  in  the 
heterozygous  variant  group,  meaning  that  we  could  not 
make  a  definitive  conclusion  regarding  the  importance 
of  heterozygous  inheritance.  By  contrast,  in  this  study, 
we  analysed  191  heterozygotes,  with  most  of  our  findings 
suggesting  that  heterozygous  inheritance  is  biologically 
relevant.  Progression-free  survival  in  the  post-prosta- 
tectomy  validation  cohort  (figure  IB)  is  the  only  set  of 
Kaplan-Meier  curves  among  the  six  plots  that  does  not 
suggest  a  difference  between  the  heterozygous  and 
homozygous  wild-type  genotypes,  and  this  cohort  had 
the  shortest  follow-up.  Apart  from  the  question  of 
whether  or  not  an  intermediate  increase  in  enzyme 
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concentration  is  associated  with  a  single  variant  allele, 
our  previous  data20  showed  that  loss  of  heterozygosity  of 
the  wild-type  allele  occurs  under  selective  pressure  from 
ADT.  Three  (27%)  of  11  patients  with  germline  hetero¬ 
zygous  inheritance  had  developed  castration-resistant 
prostate  cancer  tumours  with  loss  of  heterozygosity  of 
the  wild-type  allele,  whereas  no  patient  had  lost  the 
variant  allele.20  Loss  of  heterozygosity  of  the  wild-type 
allele  in  tumours  could  probably  contribute  to  the 
difference  in  outcomes  between  men  with  the 
heterozygous  genotype  and  those  with  the  homozygous 
wild-type  genotype.  The  discordance  between  germline 
inheritance  and  tumour  DNA  is  indicative  of  somatic 
alterations,  which  are  known  to  increase  as  tumours 
evolve,  especially  under  selective  pressure  from  ADT. 
Our  previous  work20  showed  that  even  homozygous 
wild-type  men  can  eventually  acquire  the  variant  allele 
(three  [12%]  of  25  men  assessed).  Tumour  DNA  from 
castration-resistant  prostate  cancer  is  difficult  to  obtain 
in  large  numbers  of  patients  and  was  unavailable  from 
the  cohorts  in  this  study.  Nonetheless,  we  anticipate  that 
future  studies  will  take  somatic  alterations  in  HSD3B1 
into  account  with  analysis  of  tumour  DNA  to  demonstrate 
an  even  stronger  association  between  genotype  and 
clinical  outcomes. 

Since  Huggins  and  Hodges  published  their  seminal 
work  on  the  therapeutic  effects  of  castration  in  1941, 1  ADT 
has  been  the  cornerstone  of  systemic  therapy  for  prostate 
cancer.  So  far,  prediction  of  innate  resistance  to  ADT  has 
been  challenging.  A  previous  study29  analysed  109  SNPs 
in  several  steroidogenic  genes,  including  HSD3B1,  and 
reported  no  clear  association  with  biochemical  recurrence 
in  men  with  resected  prostate  cancer.  However,  most 
patients  analysed  in  the  primary  cohort  in  that  study29 
were  not  treated  with  ADT.  This  difference  is  crucial  in 
relation  to  our  present  study,  since  the  growth  advantage 
of  the  variant  allele  would  be  most  pronounced  under 
conditions  of  castration,  as  previously  noted.  Furthermore, 
HSD3B1  has  one  homologue  (HSD3B2)  and  four  non- 
processed  pseudogenes,  which  have  very  closely  related 
DNA  sequences  that  can  obscure  detection  of  the  variant 
sequence.  The  accuracy  of  high-throughput  genotyping 
techniques,  including  the  methods  used  in  multi-SNP 
analyses,  can  be  adversely  affected  by  such  similar 
sequences  and  is  highly  dependent  on  optimal  primer 
design.30,31  Such  concerns  provided  the  motivation  for  our 
development  and  validation  of  the  specific  high-resolution 
melting  analysis  technique  we  used  in  this  study,  since 
accurate  genotyping  was  crucial. 

Another  previous  study32  associated  a  non-coding  SNP 
13  kb  upstream  of  HSD3B1  (rsl856888)  with  development 
of  castration-resistant  prostate  cancer.  This  SNP  might 
affect  the  expression  of  HSD3B1  or  have  an  independent 
regulatory  role.  Alternatively,  the  correlation  between  this 
SNP  and  castration-resistant  prostate  cancer  might  be 
due  to  its  proximity  to  HSD3B1  and  the  biochemical 
activity  that  HSD3B1  confers.  Several  other  studies33"35 


have  shown  an  association  between  germline  variants  in 
transmembrane  steroid  transporters  and  the  development 
of  castration-resistant  prostate  cancer.  Therefore,  other 
germline  variants  may  confer  additional  information  in 
combination  with  HSD3B1  (1245C). 

The  ramifications  of  a  biomarker  able  to  predict  ADT 
resistance  are  far-reaching.  Our  findings  suggest  HSD3B1 
genotype  could  a  priori  distinguish  men  with  disease  likely 
to  respond  favourably  to  ADT  from  those  with  disease 
prone  to  behave  more  aggressively,  and  who,  therefore, 
might  merit  escalated  therapy.  Combined  androgen 
blockade  with  the  upfront  addition  of  an  androgen  receptor 
antagonist  to  castration  has  long  been  debated,  but 
analyses  of  several  trials  suggest  that  the  clinical  benefit  of 
combined  androgen  blockade  in  unselected  patients  is 
small  at  best.36  Although  speculative,  a  differential  benefit 
might  be  possible  from  use  of  combined  androgen 
blockade,  with  litde  incremental  gain  for  men  who  are 
homozygous  wild-type,  but  meaningful  use  in  men 
possessing  the  variant  allele.  HSD3B1  (1245C)  genotype 
could  also  guide  future  studies  in  terms  of  selective  early 
incorporation  of  highly  potent  inhibitors  of  the  androgen 
receptor  axis,  such  as  enzalutamide  or  abiraterone  acetate. 
Men  with  the  HSD3B1  (1245C)  allele,  especially  two  copies 
(homozygotes),  could  potentially  benefit  greatly  if  one  of 
these  drugs  were  started  with  ADT  rather  than  waiting 
until  development  of  castration-resistant  prostate  cancer. 
With  respect  to  the  possible  use  of  such  drugs  for  men  who 
inherit  the  HSD3B1  (1245C)  allele,  abiraterone  is  clinically 
converted  by  3|3HSD  to  the  more  potent  anti-androgen 
metabolite  D4A,37  which  in  turn  is  further  converted  to 
other  metabolites  in  patients,  including  3-keto-5a- 
abiraterone  which  stimulates  the  androgen  receptor.38 
Nonetheless,  the  ultimate  clinical  effect  of  these  steroidal 
metabolites  of  abiraterone  remains  to  be  determined. 

Similar  reasoning  would  suggest  HSD3B1  genotype 
could  be  informative  in  decisions  regarding  chemo- 
hormonal  therapy.  Two  landmark  trials,  the  Eastern 
Cooperative  Oncology  Group  E3805  (CHAARTED)39  and 
STAMPEDE,*  showed  that  use  of  chemohormonal 
therapy  before  development  of  castration-resistant 
prostate  cancer  substantially  improved  survival  outcomes 
compared  with  ADT  alone.  Cytotoxic  therapy  would 
probably  be  most  beneficial  in  men  who  are  least  likely  to 
have  a  durable  response  to  ADT,  whereas  men  who  are 
likely  to  have  a  sustained  response  to  ADT  might  benefit 
to  a  lesser  extent.  HSD3B1  genotype  could  therefore  help 
to  guide  patient  management,  especially  for  men  with  a 
marginal  ability  to  tolerate  chemotherapy. 

The  HSD3B1  genotype  might  be  influential  in  several 
other  domains,  one  of  which  is  the  combination  of  ADT 
with  radiotherapy.  Additionally,  there  are  potential 
implications  for  more  refined  identification  of  men  for 
active  surveillance.  Active  surveillance  might  have  an 
increased  risk  in  men  with  homozygous  variant  HSD3B1 
(1245C)  inheritance,  since  they  are  not  likely  to  respond 
durably  to  ADT,  which  would  be  their  recourse  if  the 
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window  for  curative  local  therapy  were  missed. 
By  contrast,  homozygous  wild-type  inheritance  might 
reinforce  a  recommendation  for  active  surveillance. 
Similarly,  HSD3B1  genotype  could  also  be  informative  in 
consideration  of  salvage  radiotherapy  for  patients. 
Our  results  also  provide  insight  regarding  the  potential 
effect  of  pharmacological  inhibition  of  3f5HSDl.  Because 
of  the  quite  low  prevalence  of  the  homozygous  variant 
genotype,  HSD3B1  (1245C)  as  a  biomarker  could  be  used 
in  the  clinic  as  a  binary  factor — no  variant  alleles  versus 
one  or  more  variant  alleles.  We  believe  it  was  informative 
to  analyse  the  three  groups  of  genotypes  separately,  since 
the  patients  with  homozygous  variant  genotype  seem  to 
fare  especially  poorly  in  terms  of  their  rate  of  progression, 
metastasis,  and  death;  however,  we  acknowledge  that  the 
clinical  implications  might  be  similar  for  one  versus 
two  variant  alleles  in  some  clinical  situations. 

Our  study  has  several  limitations.  Each  of  the 
three  cohorts  had  a  modest  sample  size,  although  the 
similarity  of  results  in  the  three  cohorts  and  the  evidence 
of  allele  dose  dependence  strongly  support  the  effect  of 
HSD3B1  (1245C)  inheritance.  Nonetheless,  future  work 
to  prospectively  assess  the  role  of  HSD3B1  will  be  a  high 
priority,  in  view  of  the  retrospective  design  of  our  study. 
Correlation  of  serum  and  tissue  steroid  profiles  with 
genotype  will  also  be  informative.  The  prevalence  of  each 
genotype  varied  across  our  three  cohorts  (homozygous 
wild-type  37-56%,  heterozygous  variant  36-52%,  and 
homozygous  variant  6-10%).  Apart  from  random 
variation,  the  allelic  distribution  might  be  somewhat 
dilferent  in  the  populations  sampled,  even  in  the  subsets 
of  white  patients.  Future  studies  are  needed  to  provide 
additional  information  to  refine  prevalence  estimates. 
Further  analysis  of  patients  of  other  ethnic  origins  would 
also  be  of  value  in  view  of  the  under-representation  of 
these  men  in  our  cohorts,  although  the  lower  allelic 
frequency  in  non- white  populations  might  be  challenging. 

In  this  study  we  evaluated  two  cohorts  with  biochemical 
failure  and  found  HSD3B1  genotype  to  be  associated  with 
outcomes  in  this  setting.  However,  we  were  only  able  to 
meaningfully  analyse  distant  metastasis-free  survival  in 
the  primary  cohort  because  of  short  follow-up  duration  in 
the  post-prostatectomy  validation  cohort.  We  also  explored 
whether  men  at  the  opposite  end  of  the  disease  scale 
(ie,  the  validation  cohort  of  men  who  ultimately 
developed  metastatic  castration-resistant  prostate  cancer) 
had  different  outcomes  according  to  HSD3B1  status,  and 
our  data  showed  that  they  did.  However,  the  more  general 
question  of  whether  the  same  is  true  of  an  unselected 
population  of  men  with  metastatic  prostate  cancer 
warrants  further  investigation. 

One  other  potential  limitation  of  our  study  was  that  a 
uniform  follow-up  schedule  was  not  used.  Systematic 
dilferences  in  assessment  intervals  can  lead  to  spurious 
results  (eg,  in  comparison  of  progression-free  survival 
across  trials  with  different  assessment  schedules).41,42 
However,  in  this  study,  assessment  frequency  was  at  the 


discretion  of  the  treating  physicians  for  all  patients 
without  regard  to  genotype  (genotype  was  unknown  at 
the  time  of  clinical  decision  making).  Thus,  although  this 
variation  in  assessment  frequency  should  be  considered 
when  comparing  our  cohorts  with  each  other  or  with 
additional  cohorts,  it  would  not  reduce  the  validity  of 
comparing  men  by  genotype  within  each  of  our  study 
cohorts,  which  was  our  objective. 

In  conclusion,  our  findings  in  three  independent 
cohorts  nominate  germline  HSD3B1  genotype  as  a 
genetic  biomarker  of  resistance  to  ADT  for  prostate 
cancer,  and  are  concordant  with  what  would  be  predicted 
from  the  underlying  molecular  mechanism. 
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Supplementary  Methods 

DNA  extraction,  amplification,  and  genotyping. 

The  DNA  for  the  primary  cohort  was  extracted  from  3  freshly  cut  10pm  sections  of  benign  prostatectomy  tissue  using  QIAamp  DNA 
FFPE  Tissues  kit  (Qiagen,  Valencia,  CA,  USA)  following  the  manufacturer’s  instructions.  The  DNA  quantity  and  quality  was 
measured  using  NanoDrop  (Thermo  Scientific,  Waltham,  MA,  USA).  DNA  from  the  Mayo  cohorts  was  extracted  from  peripheral 
blood  mononuclear  cells.  Genomic  DNA  was  stored  at  no  more  than  4°C.  The  PCR  was  performed  in  20 llL  reactions  (final  volume) 
using  LightScanner  Master  Mix  (Biofire,  Salt  Lake  City,  UT,  USA),  4uL  extracted  genomic  DNA  ( 1  Oug/pL).  forward  primer,  reverse 
primer,  and  unlabeled  LNA  probe  at  concentrations  of  0  05,  0-5,  and  0-5pM,  respectively  (Table  SI).  The  PCR  and  the  melting 
analyses  were  performed  in  StepOnePlus™  System  (Life  Technologies,  Grand  Island,  NY,  USA)  using  the  following  conditions:  2 
min  at  95°C,  followed  by  60  cycles1,2  of  94°C  for  30  s,  66°C  for  30  s,  and  75°C  for  30  s.  The  amplification  cycles  were  then  followed 
by  the  melting  steps:  95°C  for  30  s,  followed  by  cooling  to  25°C  for  60  s,  and  a  slow  final  denaturation  to  95°C  at  a  thermal  ramp  rate 
of  0-3%.  Data  were  analyzed  using  probe  melting  derivative  peaks  (Figure  SI).  Each  validation  run  was  performed  with  duplicate 
control  samples  representing  the  three  different  genotypes  and  each  sample  was  tested  in  duplicate  during  the  study  genotyping. 
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Figures 


Melt  Curve 


Temperature  (°C) 

Figure  SI:  Derivative  melting  curves  of  homozygous  wild-type  (red),  homozygous  variant  (blue),  and 
heterozygous  (yellow).  In  a  homozygous  variant  sample,  the  probe  shows  a  melting  peak  at  a  temperature  that  is  2°C 
lower  than  the  probe’s  melting  temperature  in  a  wild-type  sample.  The  heterozygous  sample  shows  two  melting  peaks, 
whilst  the  full  length  product  in  all  samples  melts  at  a  higher  temperature  (85°C-90°C). 
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A  Distant  Metastasis-Free  Survival  According  to  HSD3B1  Genotype  for 


Mayo  SPORE  Post -Prostatectomy  Cohort 


Homozygous  Wild-Type  N=77  57  42  22  6 

Heterozygous  N=50  36  25  14  4 

Homozygous  Variant  N=10  7  2  1  1 

Overall  Survival  According  to  HSD3B1  Genotype  for  Mayo  SPORE 


Homozygous  Wild-Type  N=77  66  52  31  17  6  3 

Heterozygous  N=50  41  29  17  7  6  3 

Homozygous  Variant  N=10  7  4  1  1  0  0 


Figure  S2:  Distant  metastasis-free  survival  (A)  and  overall  survival  (B)  in  the  Mayo 
SPORE  post-prostatectomy  cohort  according  to  HSD3B1  genotype. 

Median  follow-up  in  this  cohort  (4-9  years)  was  shorter  than  for  the  other  two  cohorts,  and 
was  particularly  limited  among  the  homozygous  variant  group  (3 -2  years).  In  view  of  the 
small  number  of  events  for  DMFS  (62)  and  OS  (40)  during  the  short  observation  period,  the 
SPORE  cohort  was  used  to  validate  the  impact  of  HSD3B1  genotype  on  progression-free 
survival  only  (see  text).  Thus,  the  DMFS  and  OS  curves  above  are  presented  for  qualitative 
purposes  only,  as  the  data  are  not  mature  enough  for  reliable  quantitative  analysis  of  these 
endpoints.  The  above  limitations  notwithstanding,  the  homozygous  variant  group  does  appear 
to  diverge. 
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Progression-Free  Survival  (N=118) 


No.  at  Risk  118 


35 


2 


Distant  Metastasis-Free  Survival  (N=118) 


No.  at  Risk  118 


58  23  4 


Figure  S3:  Progression-free  survival  (A),  distant  metastasis-free 
survival  (B),  and  overall  survival  (C)  in  the  primary  study  cohort  (all 
genotypes  combined). 
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Figure  S4.  Progression-free  survival  in  (A)  post-prostatectomy  and 
(B)  metastatic  validation  cohorts  (all  genotypes  combined). 
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Overall  Survival  in  a  Metastatic  Validation  Cohort  (N=188) 


No.  at  Risk  188  169  130  93 


Figure  S5:  Overall  survival  in  a  metastatic  validation  cohort. 


Prostate  Cancer-Specific  Survival  According  to  HSD3B1  Genotype  (N=118) 


No.  at  Risk 
Homozygous  Wild-Type 
Heterozygous 
Homozygous  Variant 


44 

62 

12 


36 

53 

11 


24 

35 

7 


11 

15 

0 


Figure  S6:  Prostate  cancer-specific  survival  in  the  primary  cohort  according 
to  HSD3B1  genotype.  Prostate-cancer  specific  survival  diminished  according  to 
the  number  of  variant  HSD3B1  alleles  inherited  (P=0  029). 
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Tables 


Oligonucleotide  Sequence 

Forward  primer  5’-GT CAAAT AGCGT ATT CACCTT CT CTT AT-3’ 

Reverse  primer  5-GAGGGTGGAGCTTGATGACATCT-3’ 

Unlabeled  probe*  5’-GGAG AACCT GAAGTCCAAGACT CAGT GATTT AAGG-3’ 

*The  underlined  bold  character  indicates  the  position  of  the  mutation  and  the  LNA  base. 

Table  SI:  Oligonucleotide  sequences  used  in  this  study 


Characteristic 

Homozygous  wild-type 

N  =  77/137  (56%) 

Heterozygous 

N  =  50/137  (37%) 

Homozygous  variant 

N  =  10/137  (7%) 

P-value 

Age  (years) 

Median 

67 

66 

68 

0-47 

Interquartile  range 

61-72 

61-70 

64-75 

Race  -  no.  (%) 

White 

75  (97) 

47  (94) 

10  (100) 

0-52 

Unknown 

1  (1) 

0 

0 

Other 

1  (1) 

3(6) 

0 

T  stage  -  no.  (%) 

T2 

35  (46) 

20  (40) 

2  (20) 

0-20 

T3a 

8  (10) 

12  (24) 

1  (10) 

T3b4 

24  (31) 

10  (20) 

2(20) 

TX 

10(13) 

8(16) 

5(50) 

N  stage  -  no.  (%) 

NO 

67  (87) 

42  (84) 

5(50) 

0-025 

N1 

10(13) 

8(16) 

5  (50) 

PSA  at  ADT  initiation  (ng/mL) 

Median 

0-7 

0-7 

1-2 

Interquartile  range 

0-3-20 

0-3-2-2 

0-8-3- 1 

Gleason  score  -  no.  (%) 

<6 

5(7) 

4(8) 

0 

0-99 

7 

31  (40) 

21  (42) 

5  (50) 

8-10 

37  (48) 

23  (46) 

5  (50) 

Unknown 

4(5) 

2(4) 

0 

Neoadjuvant  ADT  -  no.  (%) 

Y 

19(25) 

6(12) 

2(20) 

0-20 

N 

58  (75) 

44  (88) 

8(80) 

Adjuvant/Salvage  RT  -  no.  (%) 

Y 

31  (40) 

29  (58) 

6(60) 

0-97 

N 

46  (60) 

21  (42) 

4(40) 

SPORE  denotes  the  National  Institutes  of  Health-funded  Specialized  Programs  of  Research  Excellence.  PSA  denotes  prostate- 
specific  antigen  level.  Staging  refers  to  the  American  Joint  Commission  on  Cancer  Tumor-Node-Metastasis  (TNM)  staging 
system.  ADT  denotes  androgen  deprivation  therapy.  P-values  represent  comparisons  across  the  three  genotypes. 


Table  S2:  Demographic  and  treatment  characteristics  of  the  Mayo  Clinic  SPORE  cohort. 
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Characteristic 

Homozygous  wild-type 

N  =  98/188  (52%) 

Heterozygous 

N  =  79/188  (42%) 

Homozygous  variant 

N  =  11/188  (6%) 

P-value 

Age  (years) 

0-62 

Median 

65 

67 

65 

Interquartile  range 

60-70 

60-73 

57-68 

Race  -  no.  (%) 

0-74 

White 

94  (96) 

76  (96) 

11  (100) 

Unknown 

1  (1) 

2(2.5) 

0 

Other 

T  stage  -  no.  (%) 

0-068 

T1-2 

43  (44) 

28  (35) 

1  (9) 

T3-4 

45  (46) 

39  (49) 

6(55) 

TX 

10  (10) 

12  (15) 

4(36) 

N  stage  -  no.  (%) 

0-38 

NO 

41  (42) 

32  (41) 

3(27) 

N1 

27  (28) 

21  (27) 

6(55) 

NX 

30 (31) 

26  (33) 

2(18) 

M  stage  at  diagnosis  -  no  (%) 

0-43 

MO 

7(7) 

11  (14) 

0 

Ml 

31  (32) 

22  (28) 

5(46) 

MX 

60  (61) 

46  (58) 

6  (55) 

PSA  at  ADT  initiation  (ng/mL) 

0-21 

Median 

7-5 

120 

39-9 

Interquartile  range 

1  -6-28-1 

2-8-61 -8 

2  0-266-3 

Gleason  score  -  no.  (%) 

0-42 

<6 

9(9) 

6(8) 

0 

7 

33  (34) 

25  (32) 

1  (9) 

8-10 

52  (53) 

41  (52) 

9  (82) 

Unknown 

4(4) 

7(9) 

1  (9) 

Chemotherapy  -  no.  (%) 

0-83 

Y 

56  (57) 

48  (61) 

6  (55) 

N 

42  (43) 

31  (39) 

5  (46) 

PSA  denotes  prostate-specific  antigen  level.  Staging  refers  to  the  American  Joint  Commission  on  Cancer  Tumor-Node- 
Metastasis  (TNM)  staging  system.  ADT  denotes  androgen  deprivation  therapy.  P-values  represent  comparisons  across  the  three 
genotypes. 


Table  S3:  Demographic  and  treatment  characteristics  of  the  Mayo  Clinic  metastatic  cohort. 
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Progression-free 

survival 

Distant  metastasis- 

free  survival 

Overall  survival 

Characteristic 

HR  (95%  Cl) 

P-value 

HR  (95%  Cl) 

P-value 

HR  (95%  Cl) 

P-value 

HSD3B1  genotype 

Homozygous  variant  vs.  homozygous  WT 

1-9  (0-8  to  4-2) 

0-13 

2-8  (1-1  to  6-7) 

0  025 

3-5  (1-3  to  9-5) 

0013 

Heterozygous  vs.  homozygous  WT 

1-6  (0-9  to  2-7) 

0-10 

1-8  (10  to  3-3) 

0  050 

2  0  (10  to  3-9) 

0  054 

At  least  one  variant  allele  vs.  none* 

1-6  (10  to  2-7) 

0  074 

T  stage 

T3a  vs.T2 

2-5  (1-1  to  5-6) 

0  033 

1-1  (0-5  to  2-5) 

0-74 

10  (0-4  to  2-5) 

0-94 

T3b-4  vs.  T2 

2-5  (1-1  to  5-9) 

0  030 

1-1  (0-5  to  2-4) 

0-77 

1-5  (0-6  to  3-4) 

0-40 

T3b-4  vs.  T3a 

10  (0-6  to  1-9) 

0  91 

10  (0-5  to  1-9) 

0-96 

1-5  (0-7  to  3-1) 

0-27 

N  stage 

N1  vs.  NO  or  NX 

1-7  (10  to  3  0) 

0  062 

1-4  (0-8  to  2-5) 

0-30 

1-5  (0-8  to  2-7) 

0-21 

PSA  at  ADT  initiation 

>2  ng/mL  vs.  <2  ng/mL 

2-4  (1-5  to  4  0) 

0  0003 

2-2  (1-3  to  3-6) 

0  0032 

1-9  (1-1  to  3-4) 

0  022 

Gleason  score 

8-10  vs.  6-7 

1-3  (0-8  to  2  0) 

0-36 

1-8  (1-1  to  3  0) 

0  023 

1-7  (0-9  to  2-9) 

0  076 

HR  denotes  hazard  ratio.  Cl  denotes  confidence  interval.  WT  denotes  HSD3B1  homozygous  wild-type  genotype.  PSA  denotes 
prostate-specific  antigen  level.  Staging  refers  to  the  American  Joint  Commission  on  Cancer  Tumor-Node-Metastasis  (TNM) 
staging  system.  ADT  denotes  androgen  deprivation  therapy.  ‘When  we  evaluated  the  impact  of  inheriting  at  least  one  variant  allele 
on  progression-free  survival,  the  adjusted  outputs  for  covariates  were  nearly  identical  and  thus  the  table  was  not  duplicated  (data 
not  shown).  ■■  =  not  assessed. 

Table  S4:  Multivariable  analyses. 
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Distribution  of  Progression  Events 

As  mentioned  in  the  text,  progression-free  survival  was  a  composite  endpoint.  For  the  primary  study  cohort,  the  distribution  of  the 
component  types  of  events  was  recorded.  Of  initial  progression  events,  70%  were  based  on  PSA  criteria  (second  increase  on  ADT); 
many  of  these  men  subsequently  progressed  clinically  (ultimately  reflected  in  the  distant-metastasis  free  survival  and  overall  survival 
curves).  The  remaining  30%  of  initial  progression  events  consisted  of  22%  that  were  based  on  clinical  progression,  and  8%  based  on 
initiation  of  second-line  therapy. 
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